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A DYNAMIC  METHOD  OF  DETERMINING  THE  STIFFNESS  AND  CROSS  AXIS 


STABILITY  OF  A REPULSION  MAGNETIC  BEARING 


SUMMARY 


Magnet  bearing  support  systems  are  becoming  of  increasing  interest  in 
satellite  and  other  engineering  projects.  In  order  to  design  such  systems  a 
knowledge  of  the  bearing  stiffness  is  required.  This  Report  analyses  the  dynamics 
of  a repulsive  type  magnet  bearing  and  proposes  a simple  pendulum  experiment  to 
determine  the  radial  stiffness  by  measuring  the  vibration  frequency.  The 
analysis,  based  on  the  potential  energy  of  the  system,  shows  the  relationship 
between  radial  stiffness  and  the  pendulum  arm  length  and  also  predicts  a rota- 
tional stability  condition,  namely  that  the  bearing  length  to  diameter  ratio  must 
be  greater  than  /3  for  stability.  ^ 
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1 INTRODUCTION 

In  recent  years  there  has  been  increasing  interest  in  the  use  of  magnetic 
bearings  for  rotating  mechanisms,  such  as  momentum  wheels  in  satellites,  in  order 
to  obtain  a reduction  in  friction  and  wear  of  conventional  bearings  and  hence  a 
greater  reliability  and  lifetime.  In  order  to  design  such  bearings  it  is  neces- 
sary to  determine  the  stiffness  of  the  system.  One  way  of  doing  this  is  simply 
to  directly  measure  the  forces  involved  with  a static  measuring  system  using 
force  transducers  and  position  measuring  equipment . Such  a procedure  introduces 
measuring  uncertainties  due  to  equipment  biases  and  tends  not  to  give  very 
accurate  results.  The  method  is  even  less  accurate  when  attempting  to  estimate 
the  cross-axis  stability  due  to  the  angular  stiffness. 

An  alternative  method  is  proposed  here,  with  the  appropriate  theoretical 
analysis,  based  on  a pendulum  type  experiment.  Not  only  does  this  method  provide 
the  stiffness  of  a bearing  but  it  can  also  be  used  to  confirm  the  'cockling  ratio' 
which  is  predicted  by  the  analysis:  namely,  that  a bearing  is  angularly  stable 
if  its  length  is  greater  than  /3  times  its  diameter.  This  has  been  confirmed  by 
the  experimental  work  of  Lodge. 

2 

The  only  theoretical  analysis  that  is  known  to  the  author  is  that  of 
Backers'.  This  work  predicted  the  positional  stiffness  of  a magnet  but  not  the 
angular  stiffness  and  was  limited  to  a two-dimensional  analysis.  This  latter 
restriction  is  removed  in  the  present  paper  where  a full  three-dimensional 
analysis  is  made. 

2 DESCRIPTION  OF  THE  EXPERIMENT 

The  type  of  repulsion  magnetic  bearing  considered  here  is  illustrated  in 
Fig  1 . Such  a bearing  consists  of  an  outer  shell  of  permanent  magnet  rings  of 
equal  thickness,  magnetised  axially  with  the  rings  arranged  so  that  like  poles 
are  adjacent,  ie  alternate  rings  have  opposite  directions  of  magnetisation. 

Inside  this  is  a coaxial  inner  shell  of  magnetic  rings  of  the  same  thickness  and 
stacked  on  a shaft  in  exact  correspondence  with  the  outer  shell.  Such  an  arrange- 
ment is  radially  stable  but  axially  unstable  and  the  angular  stability  depends, 
as  mentioned  above,  on  the  dimensions  of  the  bearing. 

If  the  shaft  of  the  bearing  is  now  suspended  from  a point  S,  Fig  2,  on  the 
bearing  axis  whilst  the  outer  shell  is  fixed  vertically  the  inner  shell  is  free 
to  oscillate  and  the  frequency  of  oscillation  will  be  a measure  of  the  forces 
and  couples  acting  on  the  bearing.  We  will  now  determine  what  this  relationship 
is. 
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3 THE  FIELD  OF  THE  OUTER  SHELL 

2 

Following  Backers  we  will  assume  that  the  intensity  of  magnetisation  M 
of  the  magnets  follows  a sinusoidal  spatial  distribution  axially,  i-e 

M - V sin  2TrZ/X  (1) 


where  X is  the  periodicity  of  the  system,  so  that  each  magnet  is  of  thickness 
t ■ X/2  , and  ^ is  the  maximum  magnetisation  occurring  at  the  centre  of  a 
magnet.  The  Z axis  is  taken  along  the  bearing  axis  with  origin  at  the  inter- 
face between  two  adjacent  magnets.  It  will  also  be  assumed  that  the  magnets  are 
homogeneous  so  that  the  field  is  radially  symmetric  and  end  effects  will  be 
ignored . 

The  magnetic  potential  may  now  be  calculated  from  the  formula 

where  the  integration  is  over  the  outer  shell  or,  integrating  by  parts. 


and  the  second  integral  is  zero  since  the  surface  element  d£  is  normal  to  the 
magnetisation  vector.  Consequently,  the  magnetic  potential  at  a point  Zq  along 
the  bearing  axis  is 


“<V 


.4'a  cos  aZ 

pdp  ■ - 

+ (Z  - Zq)^ 


(2) 


where  and  p^  denote  the  inner  and  outer  radii  respectively  of  the  outer 

shell  and  a * 2it/X  “ ir/t  . The  integral  in  equation  (2)  may  be  evaluated 
explicitly  to  give 


^(Zq)  - - 4ir.<rcos  aZQ[p,Kj  (apj^)  - PqKj  (ap^)] 


(3) 
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where  K^(x)  denotes  the  modified  Bessel  function  of  the  second  kind  defined  by 


K^(x)  ■ j e”*  ^ cosh  nt  dt  . (4) 

0 

Thus  the  magnetic  potential  along  the  bearing  axis  varies  as 

S2(Z}  - F cos  aZ  (5) 

where  F is  a negative  constant  defined  by  equation  (3)  and  the  suffix  0 has 
now  been  omitted.  The  explicit  form  for  F will  not  be  needed  in  this  Report 
since  it  is  not  necessary  to  evaluate  the  absolute  forces  and  couples  but  only 
their  relative  magnitudes. 

Now  the  magnetic  potential  must  satisfy  Laplace's  equation  and  consequently 
its  spatial  distribution  within  the  region  occupied  by  the  inner  shell  is 

■ F cos  aZ  IqCuR)  (6) 

where  R denotes  the  radial  position  from  the  bearing  axis  and  I (x)  denotes 

n 

the  Bessel  function  with  purely  imaginary  argument,  defined  by 

Iq(x)  " ^ ^ • (7) 

0 

Thus  the  magnetic  field  ^ grad  n is 


H 


- oF  sin  aZ  lQ(aR)Z  - oF  cos  aZ  Ij (aR)R 


(8) 


where  ^ and  R denote  unit  vectors  in  the  axial  and  radial  directions 
respectively.  Fig  3 shows  the  distribution  of  the  magnetic  field  and  potential. 

4 FCTICES  AND  COUPLES  ACTING  ON  A SINGLE  INNER  RING 

It  is  shown  in  the  Appendix  that  for  small  axial  and  radial  displacements 
of  the  centre  of  an  inner  ring,  a and  d respectively,  and  a small  angular 
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rotation  6 about  the  equilibrium  position  the  potential  energy  of  a single  ring 
is  approximately 


W - Wq  j^cos  aa  I^Cod)  - f (c^r^  - I-  + 

where  is  the  outer  radius  of  an  inner  ring.  Consequently  the  axial  and 

radial  forces  F and  F acting  on  the  centre  of  an  inner  ring  are 

ft  IT 

arj 

Fa  - “ 37  ’ * “^0  ““ 

JiU 

F^  - “ ad"  “ ” Ij(ad) 


which  approximate  for  small  displacements  to 


(9) 


2 

where  S.  ■ a is  the  axial  stiffness.  It  is  to  be  noted  that  this  shows  that 
A 0 

the  axial  stiffness  is  twice  the  radial  stiffness  independently  of  any  small 

rotation  present. 

Similarly,  for  small  displacements,  the  couple  acting  about  the  centre  of 
a ring  is 


C 


IE 

3e 


(10) 


It  is  to  be  noted  that  this  couple  has  a destabilising  influence. 

5 MOTION  OF  THE  PENDULUM 

Let  us  now  consider  the  total  couple  acting  on  the  inner  shell  assembly, 
about  the  point  of  suspension.  If  there  are  a total  of  n inner  magnets  so  that 
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the  length  of  the  inner  assembly  is  » nt  and  L is  the  distance  from  the 
point  of  suspension  to  the  centre  of  the  inner  assembly  then  the  total  moment 
about  the  point  of  suspension  due  to  the  radial  forces  acting  on  each  magnet  is, 
using  equation  (10) 


“ - S„e  I [l  + (N  + i - n)t]' 

r tv  1 


-S^en 


-zJ-t^' 
12  ^ 


The  moment  due  to  the  axial  force  is  second  order  in  6 since  the  axial  displace- 
ment depends  upon  1 - cos  6 and  may  thus  be  neglected  to  first  order  in  6 . 
Adding  now  the  destabilising  couple  of  each  magnet  gives  a total  moment  about 
the  point  of  suspension 


M + nC 
r 


-S^en 


where  D - 2rQ  is  the  diameter  of  an  inner  ring. 

From  this  equation  it  can  be  seen  that  if  there  is  a long  suspension  arm, 

then  the  moment  about  the  point  of  suspension  is  equal  to  the  arm  length  L 

multiplied  by  the  total  radial  stiffness  ■ nS^^  of  the  bearing,  ie  the 

destabling  couples  may  be  ignored.  On  the  other  hand,  when  the  inner  bearing  is 

supported  about  its  centre  (or  freely  suspended  with  the  aid  of  a stabilising  axial 

2 2 2 2 

servo),  the  moment  about  its  centre  depends  upon  the  quantity  i - 3D  - (3t  /it  ). 
2 2 

Since  the  term  3t  /n  is  negligible,  the  bearing  will  thus  be  rotationally  stable 
if  its  length  is  greater  than  ^ times  its  diameter  and  unstable  otherwise. 

If,  now,  m is  the  mass  of  the  inner  bearing  and  k is  its  radius  of 
gyration  about  an  axis  through  its  centre  normal  to  the  bearing  axis,  then  the 
total  moment  of  inertia  about  the  point  of  suspension  is 


m(L^  + k^)  + I 
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where  is  the  moment  of  inertia  of  the  supporting  pendulum  shaft.  Also, 

if  ~0Mg  is  the  moment  about  the  point  of  suspension  due  to  gravitational  forces, 
the  total  moment  is  ~ 0Mg  and  consequently  the  oscillation  frequency  of  the 

pendulim,  u , in  rad/s  is  given  by 


t2  1 

/«2  ..  2 3t^\ 

^RT 

L 

U -3D  j 

+ M 

g 

m(L^  + k^)  + I 

arm 


(13) 


The  gravitational  term  will  be  positive  if  the  point  of  suspension  is  above 
the  pendulum,  otherwise  negative,  and  is  proportional  to  the  arm  length  L . For 
practical  bearings  this  term  will  be  insignificant  compared  with  that  contributed 
from  the  radial  stiffness.  For  long  arm  lengths,  then,  the  oscillation  frequency 
will  approach  the  asymptotic  value 


CO 


a 


(14) 


independently  of  the  arm  length.  Consequently  the  period  of  oscillation  provides 
a measure  of  the  radial  stiffness  of  the  bearing  and,  by  varying  the  arm  length, 
the  rotational  stability  condition  can  be  investigated.  Fig  5 shows  the  predicted 
variation  of  co^  with  . 

6 CONCLUSION 


The  analysis  has  shown  the  relationship  between  the  radial  stiffness,  the 
geometry  of  the  bearing  suspension  system  and  the  frequency  of  oscillation, 
equation  (13).  By  measuring  the  period  of  oscillation  for  long  arm  lengths  a 
direct  measure  of  the  radial  stiffness  is  available.  Furthermore  by  varying  the 
arm  length  the  rotational  stability  condition  can  be  investigated.  A limited 
number  of  experiments  have  been  reported  in  Ref  1 and  support  the  predictions  of 
this  Report. 
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Appendix 

POTENTIAL  ENERGY  OF  A SINGLE  RING 
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Consider  the  situation  where  the  inner  shell  is  rotated  by  a small  angle 
6 about  the  point  of  suspension  S,  Fig  4.  Let  this  be  in  the  plane  XZ  with 
the  Y axis  forming  a right  hand  set,  and  let  a and  d denote  corresponding 
axial  and  radial  displacements  of  the  Nth  inner  magnet.  We  will  now  introduce 
another  cartesian  frame  of  reference  (xyz)  with  its  origin  at  the  centre  of  the 
magnet,  z axis  along  the  axis  of  the  magnet  and  the  xz  plane  coinciding  with 
the  XZ  plane.  Thus  if  (x,y,z)  are  the  coordinates  of  a point  in  the  magnet, 
then  the  corresponding  (X,Y,Z)  coordinates  are 


X ■ R cos  $ = d + X cos  0 + z sin  0 

Y = R sin  = y 

Z « (N-j[)t  + a-  x sin  0 + z cos  0 (A-1) 


where  (R,t,Z)  are  the  corresponding  cylindrical  coordinates. 

Now  the  potential  energy  W of  the  ring  is  given  by  the  equation 

W - - WqJ  H . Mj^dv  (A-2) 

where  H is  the  magnetic  field  of  the  outer  shell,  ^ is  the  magnetisation  of 
the  inner  ring  and  the  integration  is  over  the  volume  of  the  inner  ring.  Now 

M^(z)  « sin  [a(N  - i)t  + zj  £ 

■ .^cos  az  £ (A-3) 

where  £ is  a unit  vector  along  the  z axis.  Consequently,  by  integrating 
equation  (A-2)  by  parts  and  noting  that  G •»  0 at  the  z boundaries,  the  poten- 
tial energy  is 


f 3M, 


W - 
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Appendix 


f 


L 


and  thus,  using  equations  (A-1),  (A-3)  and  (6), 

~ 1^0  JdvaF^sin  oz  sin  a (a  - x sin  6 z cos  e)lQ(aR) 

Introducing  polar  coordinates  (r,4/)  , so  that  x - r cos  , y - r sin  t/;  , 

and  the  integral  representation  (7)  for  Iq  gives 

W - - jrdrdilidz  sin  az  sin  o (a  - r cos  ij;  sin  9 + z cos  0) 

2ir 

* i dt  . (A-4) 

0 


For  convenience  we  will  now  use  the  thickness  of  a magnet  as  the  unit  of  length 

so  that  we  may  put  a-1.  Using  equation  (A-1)  equation  (A-4)  may  then  be  put 
in  the  form 


Jo  ,r/2 

W - - UqF.  J'  j rdv  J dz  sin  z Im  J 


(A-5) 


r.  -7r/2 


where  Im  J denotes  the  imaginary  part  of  the  integral 


27t  2it 


■z  cos  6 - r cos  \p  sin  6) 


■ 27  J I 

X + r cos  <li  cos  6 + z sin  9)  cos  t + r sin  \|)  sin  t 

and  rQ  and  r^  denote  the  inner  and  outer  radii  of  the  inner  shell. 
Expanding  J to  second  order  in  0 about  9-0  produces 
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1 1 


2Tr  2Tr 


2'nJ 


/"  j , i(a+z)  + d cos  t + r 
■ dt  dilie 

J J 


cos 


[i  - e(i 
e2  /. 


ir  cos  ip  - 2 cos  t) 


. 2 2 2 2 ' 
12  + r cos  ill  cos  t + 2ir2  cos  ip  cos  t + r cos  ip  - 2 cos  t 


‘)] 


and  using  equation  (7)  to  evaluate  the  integrals  results  in 

i [lo(r)I;,(d)  - ei,(d)^ltl|(r)  - zl|,(d)^ 

! 'n 

- — I i2lQ(r)lQ(d)  + r(l  + 2i2)Ij(r)  — 


I^(d)+l2(d) 


■)]• 


Now  consider  the  2 integration  occurring  in  equation  (A-5)  for  W. 
Firstly,  the  term  independent  of  rotation  in  equation  (A-6)  contributes  a 
potential  energy 


Wj  = - 2iruQF.^i 


yj  II  / 

/, 


r . 

I 


7t/2 

i 


dz  sin  z sin  (a  + 2)lQ(r)lQ(d) 


- Wq  cos  a lQ(d)  (A-7) 

where  is  the  equilibrium  potential  energy  of  a ring 

“0  ■ - • <*-®> 

The  remaining  terms  in  equation  (A-6)  are  already  of  second  order  for  d small 
and  hence,  to  second  order,  we  may  put  a zero,  ie  the  rotational  energy  is 
independent  of  small  axial  motion.  Furthermore,  with  a - 0 , the  0 coeffi- 
cient occurring  in  equation  (A-6)  gives  zero  contribution  to  W since  the 
integrand  of  the  2 integration  in  equation  (A-5)  is  an  odd  function  of  z. 

From  this  we  may  conclude  that  for  small  translation  and  rotation  about  the 
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Appendix 


equilibrium  position,  the  nett  forces  acting  at  the  centre  of  an  inner  magnet  are 
independent  of  rotation  and  the  couple  is  independent  of  translation. 

2 

The  remaining  9 coefficient  in  equation  (A-6)  contributes,  for  a and 
d small,  a rotational  energy 


*^0  ir/2 


W, 


r.  -i/2  ' 


+ 2z  cos  z (ln(r)  + rlj(r))|  . 


>)j  • (A-9) 


Using  the  basic  integrals 


Jtljdr  - 

rl]  , 

tt/2 

Tr/2 

1 sin^x 

dx  - J , 

h 

sin  2x  dx 

-1/2 

-i/2 

reduces  equation  (A-9)  to 


r . 
1 


Now  provided  that  (r^  - r^^)  is  several  times  the  thickness  of  a magnet, 
as  is  always  the  case  in  practice,  the  contribution  from  r^  in  equation  (A-10) 
is  negligible.  Consequently,  the  rotational  energy  may  be  approximated  by 


e2  2 


using  equation  (A-8) . 
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Summarising,  then,  and  re-introducing  the  factor  a the  potential  energy 
of  a displaced  inner  ring  is  given  approximately  by 

W - Wq  I^cos  aa  I^Cad)  - ^a^r^  " TI  * I 

for  small  displacements  a , d and  a small  rotation  6 . 
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